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We study theoretically the properties of buckled graphene-like materials, such as silicene and 
germanene, in a strong perpendicular magnetic field and a periodic potential. We analyze how 
the spin-orbit interaction and the perpendicular electric field influences the energy spectra of these 
systems. When the magnetic flux through a unit cell of the periodic potential measured in magnetic 
flux quantum is a rational number, a = p/q, then in each Landau level the energy spectra have 
a band structure, which is characterized by the corresponding gaps. We study the dependence of 
those gaps on the parameters of the buckled graphene-like materials. Although some gaps have weak 
dependence on the magnitude of the spin-orbit coupling and the external electric field, there are 
gaps that show strong nonomonotic dependence on these parameters. For a = 1/2, the spin-orbit 
interaction also opens up a gap at one of the Landau levels. The magnitude of the gap increases 
with spin-orbit coupling and decreases with the applied electric field. 


I. INTRODUCTION 

Electron dynamics in a periodic potential subjected 
to a perpendicular magnetic field was the subject of in¬ 
tensive theoretical and experimental research for several 
decades That interest on this particular subject 

derives from its unique energy spectrum, which as a func¬ 
tion of the magnetic flux per lattice cell of the periodic 
potential reveals a fractal pattern, known as the Hof- 
stadter’s butterfly (due to the pattern resembling the 
butterflies). This is the first example of fractal pattern 
realized in a quantum system. The parameter which de¬ 
termines the fractal structure of the energy spectrum is 
the magnetic flux a through a unit cell measured in units 
of the magnetic flux quantum. In the case of a weak 
magnetic field B, first the periodic potential forms the 
Bloch bands and then the external magnetic field splits 
each Bloch band of the periodic potential into minibands 
of Landau level (LL) type, the number of which is de¬ 
termined by the parameter a. In the strong magnetic 
field regime, the relevant parameter is a = 1/d which 
is zero when B —> oo. In this case, the energy spec¬ 
tra of an electron can be described by the formation of 
the LL spectrum and then splitting of the LL states into 
minibands by the periodic potential. Here the number 
of minibands is determined by the parameter a, i.e., for 
rational a = p/q where p and q are integers, within a 
single LL there are q minibands with p degeneracy. To 
observe the fractal pattern in a reasonable range of the 
magnetic field, the important requirement is that the lat¬ 
tice structure has a large period. This can be achieved 
in artificial superlattices based on semiconductor nanos¬ 
tructures Observation of a clear fractal pattern re¬ 
mained a challenge however. Finally, in 2013 the Hofs- 
tadter butterfly structure with clear fractal pattern was 
reported [kJ-QJ] hi monolayer and bilayer graphene fl3| . 
In these experiments, the periodic lattice with a period 
ss 10 nm was created naturally by the moire pattern that 


appears when graphene is placed on hexagonal boron ni¬ 
tride with a twist 0, H8j . In addition to the natural 
formation of the period structure in graphene on boron 
nitride substrate, that system has also the unique rela¬ 
tivistic energy disp ersion, which results in special type of 
Landau levels |l4Hl6j . 

Recently discovered new Dirac-type materials, such as 
silicene and germanene [ToI - IttI bring in additional re¬ 
markable features and additional control parameters into 
the structure of the Dirac energy spectrum. These ma¬ 
terials are monolayers of silicon and germanium with 
hexagonal lattice structures where the low energy charge 
carriers are also massless Dirac fermions [28140 just as in 
graphene 0 , 0 - Experimentally, the two-dimensional 
(2D) silicene was synthesized on Ag(lll) and zir¬ 

conium diboride substrates 34], while germanene was 
grown on Ag (H, 0 and Pt [36j substrates. The main 
difference between silicene/germanene and graphene is 
that due to the larger radius of the Si/Ge atom com¬ 
pared to the C atom, the corresponding hexagon lat¬ 
tices in germanene and silicene have buckled structure 
(37j, i.e., the two sublattices (say A and B) in these sys¬ 
tems are displaced vertically by a finite distance L z . As 
a result, silicene and germanene have large spin-orbit 
interactions, which opens the band gaps at the Dirac 
points (A so ss 1.55 — 7.9 meV for silicene 0, 0 and 
A so « 24 — 93 meV for germanene 0, 0). In the case 
of graphene, in contrast, the corresponding spin-orbit- 
induced gap is tiny, 25 peV 0]. The buckled structure 
of silicene/germanene lattice also allows for the band gap 
to be controlled by an applied perpendicular electric field 
0 and the size of the band gap increases almost linearly 
with the electric field. These properties have important 
implications for a highly correlated electron systems in 
these graphene-like but novel systems [ill ], in particular, 
in the fractional quantum Hall effect regime 16,142]. The 
buckled graphene-like materials with a strong spin-orbit 
interaction and sensitivity to the external electric field 
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FIG. 1: Energy spectrum (Hofstadter butterfly) of germanene 
monolayer as a function of the parameter a , which is the 
inverse magnetic flux through the unit cell in units of the flux 
quantum. The external electric field is E z = 50 mV/A. The 
amplitude of the periodic potential is V 0 = 40 meV and its 
period is a = 20 nm. 


significantly modify the energy spectrum of the system, 
which should be also visible in the Hofstadter’s butterfly 
pattern of these materials. Below we study the proper¬ 
ties of the buckled graphene-like materials, silicene and 
germanene, placed in a perpendicular magnetic field and 
a periodic potential. 


II. THEORETICAL FRAMEWORK 


We consider an electron in a silicene/germanene mono- 
layer in an external perpendicular magnetic field and a 
periodic potential. The Hamiltonian, 'H, of such an elec¬ 
tron consists of the kinetic energy term TL 0 and the peri¬ 
odic potential V(x,y), 

H = U 0 + V(x,y). (1) 


In a magnetic field B , the kinetic energy part of the 
Hamiltonian has the following matrix form [43] 
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where A R is the Rashba spin-orbit constant, £ 0 = 
yJhc/eB is the magnetic length, w c = \/2hv F /£ 0 , v F 
is the Fermi velocity, A ± = tA so + L Z E Z . Here A so 
is the spin-orbit constant, 2 L z is the separation of two 
sublattices A and B in the z direction, and E z is the 
external perpendicular electric field. For germanene and 
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FIG. 2: Energy spectra of germanene monolayer for a = 1/3. 
The external electric field is E z = 10 mV/A. The numbers 
near the lines are the amplitudes of the periodic potential, 
Vo- The period of the potential is a = 20 nm. Two gaps A 1 
and A 2 are marked in the figure. 


silicene, the parameters in the above Hamiltonian are 
u F = 7.26 x 10 5 m/s, L z = 0.33 A, A so = 43 meV, 
Ar, = 10-7 meV for germanene and v F = 8.47 x 10 5 m/s, 
L z = 0.23 A, A so = 3.9 meV, A R = 0.7 meV for silicene. 

The wavefunctions of the Hamiltonian m have four 
components and can be expressed in terms of the nonrel- 
ativistic wave functions <f> n k , which correspond to n-th 
conventional nonrelativistic Landau level and have the 
in-plane y component of the wave vector k , 

I nk) = [C l (l) nk: C2(l) n+ i^kTC 3 (j) n _ l k ,C i (j) nk ] . (3) 

In this basis the Hamiltonian (0 takes the form 
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where K n = hu> c y/n, a n = y/2nha\ R /£ 0 and the 
wave functions are determined now by the coefficients 

(C lt C 2 ,C 3 ,C A ). 


The conventional Landau wave functions <f) n k have the 
form 


e iky e -(x-x k ) 2 /2 

(b„ k (x, y) = —■= — = H ri (x — xu), (5) 

, ’ VL y / 7r 1 / 2 £ 0 2 ra n! 

where L is the length of a sample in the y direction, k is 
the y component of the electron wave vector, x k = Mq, 
and H n [x) are the Hermite polynomials. 

Depending on the value of n, the energy spectrum of 
the Hamiltonian (|2]) has the following properties: (i) for 
n = — 1 there is only one Landau level, which is charac¬ 
terized by the wavefunction 


1-1, *) = (Mo,fc,0,0) 


( 6 ) 
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FIG. 3: Energy spectra of buckled graphene-like materials 
for q = 1/3. The external electric field is E z = 100 mV/A 
(black lines) and E z = 0 (red lines). The numbers near the 
lines are the SO interaction constants, Ago- Th e period of the 
potential is a = 20 nm and its amplitude is Vq = 20 meV. 


and has the energy of E n __ 1 = —A + (E Z ); (ii) for n = 0 
there are three different Landau levels, 

|0, k) = 0) C^o.fc) (7) 

and their energies are determined by the following cubic 
equation 


E 3 + A _E 2 - (A* + Xl + kI) E 
+ (A?A + - A^A_ - k?A_) = 0; (8) 

In the absence of the external electric field ( E z = 0), one 
of the solutions of Eq. © is E 0 = A + = — A so . The 
corresponding wavefunction is 

|0, k) = (*oq</> 0jfc , 0,0, k^o)^) (9) 

and (iii) for n > 0 there are four differ¬ 

ent Landau levels with the general structure 

( C lKk’ C 2&n+l,k’ C 3&„-l,k> C 4&n,k)- In what fol - 
lows, we mainly consider the properties of the n = — 1 
and n = 0 Landau levels. 

The periodic potential, V(x , y), in the Hamiltonian (JT]) 
is characterized by its period a 0 and amplitude V 0 . We 
assume that the potential has the following profile 
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FIG. 4: Gaps A x and A 2 , defined in Fig. [2] as functions of the 
SO interaction, A so - The period of the potential is 20 nm. 
(a) The electric field is zero and the amplitude of periodic 
potential is 10 meV. (b) The electric field is 100 mV/A and 
the amplitude of the periodic potential is 10 meV. (c) The 
electric field is 100 mV /A and the amplitude of the periodic 
potential is 20 meV. 


have the following form 
(n'k'\cos(q 0 y) \nk) = 


{^k',k+q 0 + ( _ l) n n <V,/s-g 0 } 
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( 11 ) 


V ( x, y ) = V 0 [cos(g x ®) + cos(y y y)] , (10) 

where q x = q y = q 0 = 2ir/a 0 . The periodic potential 
mixes the electron states k within a single LL, i.e., 
states with the same value of the LL index n and different 
values of k, and also mixes the states of different LLs. 
The strength of this mixing is determined by the matrix 
elements of the periodic potential V(x,y) between the 
LL states. 

The matrix elements of the periodic potential in the 
basis of the LL wave functions of the buckled materials 
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FIG. 5: Energy spectrum of the germanene monolayer as a 
function of the effective wave vector k. The external electric 
field is zero. The amplitude of periodic potential is 40 meV 
and its period is 20 nm. 
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and Q = gg^o/2, m = min(n / ,n), M = max(n / ,n). 

We study the regime of the strong magnetic field and 
consider the basis of four LLs, which correspond to 
n = — 1 and n = 0. Within these basis we construct 
the Hamiltonian matrix taking into account the matrix 
elements m and m of the periodic potential, and eval¬ 
uate the corresponding energy spectrum as a function of 
parameter a = Bag/'bg, where >!> 0 = h/{2e) is the mag¬ 
netic flux quantum. 

III. NUMERICAL RESULTS AND DISCUSSION 
A. The butterfly Energy Spectrum 

The energy spectra of buckled graphene-like materials 
in the regime of the strong magnetic field, i.e., weak peri¬ 
odic potential, is shown in Fig. [T] as a function of the pa¬ 
rameter a , which changes from 0 to 1 . We consider only 
four LLs with the label scheme shown in the figure, i.e., 
LL_ 1 x corresponds to the LL with index n = — 1 (there 
is only one LL with that index), while LL 0 1 , LL 0 2 , and 
LL 0 3 correspond to three LLs with index n = 0. The 
results in Fig. [T| correspond to germanene, which has the 
largest SO interaction constant, A so = 43 meV. The re¬ 
sults clearly illustrate the formation of the Hofstadter 
butterfly structures in each LL. For small a , there is an 
overlap of the butterfly structures of the Landau levels 
LL 0 1 and LL_ l x . Although the overlap is large, the 
coupling of LL 0 x and LL_ 1 x states is weak. This is due 
to the fact that the overlap of the corresponding wave 
functions is small. For the zero external electric field, 
this overlap is exactly zero. For a > 0.7, i.e., in a weak 
magnetic field there is an overlap of the LL 0 x and LL 0 2 
Landau levels. In this case the overlap results in a rel¬ 
atively strong coupling of the corresponding states and 
modification of the energy spectrum. 


FIG. 6: Energy gap A between two bands in LL 0 2 - (a) The 
period of the potential is 20 nm and its amplitude is 40 meV. 
The SO coupling is A so = 43 meV. (b) The electric field is 10 
mV/A. The period of the potential is 20 nm and its amplitude 
is 40 meV. 

The results shown in Fig. [T| illustrate the typical struc¬ 
ture of the energy spectra of buckled graphene-like mate¬ 
rials. Variation of the external electric field, E z . and the 
strength of the SO interaction, A so , change the energy 
spectra of the system. These changes can be described 
by studying the values of the energy gaps for rational 
a = p/q , where q — 1 gaps exist in each LL. For example, 
for a = 1/3 there are three bands and correspondingly 
two gaps in each LL, while for a = 1/2 there are two 
bands and one gap in each LL. In the case of overlapping 
LLs the gaps can also be closed. Below we study the 
properties of the energy spectra of buckled graphene-like 
materials for a = 1/3 and a = 1/2. 


B. a = 1/3 

For a = 1/3 there are three bands in each LL with 
the corresponding two gaps. In Fig.[2]the energy spectra 
for a = 1/3 are shown for four LLs and different ampli¬ 
tudes of the periodic potential, V 0 . In each LL the three 
bands are clearly visible. With increasing V 0 the widths 
of the bands and correspondingly the band gaps increase. 
For large V 0 the coupling of the states of different LLs 
modifies the band structure of the LLs. We characterize 
this effect by studying the magnitudes of the two gaps 
marked in Fig. [2] by A-,, which is the higher energy gap 
in the Landau level LL_ l 1 , and A 2 , which is the lower 
energy gap in the Landau level LL 0 . 

The effect of the SO interaction strength on the ar- 
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rangement of the LL bands is illustrated in Fig. [3] For 
small Ag 0 (< 20 meV) the LLs LL_ X 1 and LL 0 . which 
are broadened by the periodic potential, overlap that can 
close some gaps within these LLs. This corresponds to 
the case of silicene, for which A so = 3.9 meV. With in¬ 
creasing Ago, the LLs LL_ 11 and LL 0 1 become well 
separated. Compared to other LLs, the energies of the 
LLs LL_ 1X and LL 0 1 have strong dependence on the SO 
coupling. This strong dependence results in the overlap 
of the LLs LL 0 1 and LL 0 3 for a large spin orbit cou¬ 
pling and large external electric field. The overlap of the 
LLs LL 0 i and LL 0 3 is visible in Fig. [3] for A so = 100 
meV and the electric field E z = 100 mV/A, while for 
E z = 0 there is no overlap of the LLs even for large A go ■ 
In Fig. 31 the dependence of the energy gaps A 3 and A 2 , 
which are defined in Fig. [2] on the SO coupling, A so . For 
zero electric field [Fig. [3](a)] the gaps A 3 and A 2 are ex¬ 
actly the same. In this case the dependence of A 1 on the 
SO interaction strength is relatively weak. The gap A 1 
changes by only ss 15% when the coupling A so increases 
from 3 meV (silicene) to 100 meV. 

The dependence of the gaps A 3 and A 2 on A so be¬ 
comes more pronounced for large electric fields. The gaps 
A-j and A 2 in this case are not equal. The gap A 2 still has 
a very weak dependence on A so . This gap corresponds 
to the high energy gap for LL_ 1 1 . At the same time, the 
gap A 3 has a very strong dependence on A so and with 
increasing A so it is strongly suppressed [Fig. HKb)] . The 
gap A-, changes from 2 meV for small A so ~ 3 meV to 
0.8 meV for large A so ~ 100 meV. This suppression is 
due to the fact that with increasing A so the separation 
between the LLs LL 0 3 and LL 0 3 decreases (see Fig. E), 
which results in stronger coupling of the states of these 
LLs, and correspondingly a strong change in the value 
of the gap A- t . The gap, A 1; has also a nonmonotonic 
dependence on A so , which is more pronounced for larger 
amplitudes of the periodic potential, V 0 . As an exam¬ 
ple, for V 0 = 20 meV (see Fig. (He)), the gap A 3 first 
increases from small value of ~ 2.5 meV for small A so , 
then reaches its maximum of ~ 6.5 meV for A go ~ 20 
meV, and finally decreases for large A so . 

C. a = 1/2 

For a = 1/2 there are two bands in each LL. Specific to 
this case is the fact that the gap between these two bands 
in each LL is zero. A gap can however be opened by the 
Coulomb interaction between electrons [44], Ea] • In the 
case of buckled graphene-like materials there are other 
parameters, the SO coupling and the external electric 
field, that can modify the band structure of the LLs. To 
characterize the band structure of the LLs, we introduce 
an “effective 2D wave vector”, {k x , K y ), which reflects the 
periodicity of the wave functions in the reciprocal space. 
This periodicity follows from the expressions of the ma¬ 
trix elements of the periodic potential [Eqs. (fTITl - lfl^l) ]. 
The effective wave vector ( k x , K y ) has the units of length 


and are defined within an effective Brillouin zone. In Fig. 
[5] the typical energy spectrum of the germanene layer for 
a = 1/2 is shown as a function of k x and different values 
of n y . The spectrum clearly shows two bands in each 
LL. The gaps between the bands are zero for the LLs 
LL_ 11 , LLq i, and LL 0 3 , while there is a finite gap A 
for LL 0 2 ■ In Fig. [5] the gap A is shown as a function of 
the electric field and the SO coupling. With increasing 
electric field the gap decreases and finally disappears for 
large fields. As a function of the SO interaction the gap 
increases with A go , which illustrates the fact that the 
gap is due to the strong SO coupling in the system. For 
small SO coupling, A < 20 meV, the gap is zero and it 
monotonically increases with A go reaching the value of 
2.5 meV for A go = 100 meV. Therefore, opening of the 
gap for A = 1/2 can be observed only in the germanene 
monolayer which has a large SO coupling, A ~ 40 meV, 
while in the silicene monolayer with small SO coupling, 
A go ~ 3.4 meV, the gap is zero. 


IV. CONCLUSION 


The graphene-like materials, e.g., silicene and ger¬ 
manene, have strong SO interaction and the buckled 
structure that results in these materials being sensitive 
to a perpendicular electric field. In this case the main 
parameters, which determines the unique properties of 
these materials, are the SO coupling, A go , and the ex¬ 
ternal perpendicular electric field, E z . The properties of 
buckled graphene-like materials in a magnetic field and 
the periodic potential also depend on the values of these 
parameters. That dependence can be described in terms 
of the dependence of the band structure of the LLs for 
rational values of a. One of the characteristics of the 
band structure is the gap between the intra-Landau level 
bands. For a = 1/3, some gaps show strong dependence 
on both the electric field and the SO coupling. For a 
large electric field and large amplitude of the periodic 
potential, the dependence of the gap on the SO coupling 
is highly nonmonotonic. For a = 1/2, without the SO 
coupling and for all values of the electric field, all gaps 
are closed. For a large SO interaction, A go > 20 meV, 
which is realized in germanene, the gap in one of the 
LLs opens. The magnitude of the gap increases with 
A so and decreases with the electric field. Experimental 
confirmation of these predictions will provide a rare peek 
into the electronic properties of these unique materials 
with emerging properties. Finally, the influence of the 
electron-electron interaction on the butterfly gap struc¬ 
ture is an important direction for exploration in these 
and other similar [4^ Dirac materials, that has already 
seen some progress, theoretically and experimentally, in 
graphene 111 1 ( 1 . 
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